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SIMIARY 

A con 5 )rehensive aerodynamic analysis program based on linearized 
potential theory is. described. The solution treats thickness and attitude 
problems at subsonic and supersonic speeds. Three dimensional configura- 
tions with or without jet flaps having multiple non-planar surfaces of 
arbitrary planform and open or closed slender bodies of non-circular 
contour may be analyzed. Longitudinal and lateral -directional static 
and rotary derivative solutions may be generated. 

The analysis has been implemented on a time sharing system in 
conjunction with an input tablet digitizer and an interactive graphics 
input/output display and editing teiminal to maximize its responsiveness 
to the preliminary analysis problem. Nominal case computation time of 1+5 
CPU seconds on the CDC 175 for a 200 panel simulation indicates the program 
provides an efficient analysis for systematically performing various 
aerodynamic configuration tradeoff and evaluation studies. 
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INTRODUCTION 


Con^uterization of aerodynamic theory has developed to a point where 
the analysis of complete aircraft configurations by a single program is now 
possible. Programs designed for this purpose in fact currently exist, but 
are limited in scope and abound with subtleties requiring the user to be 
highly e^erienced. Many of the difficulties are attributable to the level 
of precision of the underlying theory and the. numerical sensitivity of 
the associated solution. In preliminary design stages, it is often desirable 
to accept some degree of approximation in the interest of modest tum-around 
time, reduced con^jutational costs, simplification of input, and stability 
and generality of results . The in^ortance of short elapsed time stems from 
the necessity to survey systematically a large number of candidate advanced 
configurations or major coii5)onent geometric parameters for a set of overall 
system requironents in a timely manner. Modest computational cost allows a 
greater number of configurations and/or conditions to be economically 
investigated. 

One approach in this spirit is to replace the commonly employed exact 
superposition method, \diich panels the entire aircraft surface, with 
approximations involving linearized boundary conditions and solutions of a 
local two-dimensional potential equation. In the exact theories, a 
detemination of the singularity strengths required to satisfy boundary 
conditions leads to the necessity of inverting very large matrices . The 
nature of the approximate linearized theories on the other hand substantially 
reduces the number of simultaneous equations encountered and consequently 
places far less demand iqjon computer capabilities. 

Linearized theory when combined with realistic assessments of 
limitations and estimated viscous characteristics provides a valuable tool 
for analyzing general aircraft configurations and aerodynamic interactions 
at modest attitudes for both subsonic and supersonic speeds. 
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LIST OF SYMBOLS 


A 



Ai 

flij 


b 

c 

c 

CAVG 

Cd 

Cd 

CDp 

CDmin 

Cf 


Ci 

C£ 

C£,Qn,Cn 

CL 

Cn 

Cp 

^PNET 

Cy 


Projected oblique cross section area 

Influence coefficient. Normalwash at control point i due 
to vortex panel j of unit strength 

Area of quadrilateral panel i 

Coefficients in the set of linear equations for the vortex 
panel strengths 

Reference span 

Local chord 

Reference chord 

Average chord 

Section drag coefficient 

Drag coefficient 

Flat plate skin friction drag coefficient 

Minimum drag due to lift 

Flat plate skin friction drag coefficient 

BoundarY condition for control point i 

Section lift coefficient 

Roll, pitch and yaw moment coefficients 

Lift coefficient 

Section normal coefficient 

Pressure coefficient (P-Pdo)/q 

Net pressure coefficient CP£"Pu)/ci and vortex panel strength 
Side force coefficient 
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Cy 

c* 

D 

F 
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h 

i,j,k 

K 

I 

JlCi,n) 

L 

L/d 

m 

M 

n 

p»q,r 
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LIST OF SYflBOLS (CONTINUED) 

Section momentum coefficient mV/qc 

Uco 

Drag 

Jet reaction force 

Force components for body of unit length 

Axisymmetric outer solution to potential equation 

Radius of curvature of cross sectional boundary 

Unit vectors in x,y,z direction respectively 

Drag due to lift factor or skin friction thickness correction 
factor 

Equivalent distributed sand gain height 
Effective length 

Length of segment i, i+1 of contour Cn 
Equivalent body length or geometric length 
Bocfy fineness ratio 
Mass rate of flow 
Mach number 

Moment con^Donents for body of unit length 
Unit normal 

Rolling, pitching and yawing velocity about x, y and z 
Nondimens ional angular velocities pb/2U, qc/2U and rb/2U 
Pressure 
Prandtl number 
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LIST OF SYMBOLS CCONTINUED) 


q 

r 

R 

a 

s 

S 

Sref 

T 


t/c 


u,v,w 


U 

V 

W 

x,y,z 

x,r,e 

Z 

a 


OLi 
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Free stream dynamic pressure l/2o\j2 
Recovery factor 

Unit Reynolds number or radius of curvature 
Reynolds nimiber based on [ ] 

Gas constant 
Segment arc length 

Body cross sectional area or surface area 
Reference area 

Static tenperature -°R or tangent of quadrilateral panel 
■ leading edge sweep 

Airfoil thickness ratio 

x,y,z nondimens ional conponents of perturbation velocity 
Freestream velocity 
Jet velocity 

Conplex. potential function 
Body axis coordinate system 
Cylindrical coordinate system 
Conplex number y+iz 
Angle of attack 

Local angle of attack at surface control point i 
Angle of sideslip or J\i^ 
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LIST OF SYMBOLS (CONTINUED) 


Y 

r 

6 

'ij 

6jT 

(Sv 

A(j) 

e 

u 

V 
p 

<T 

♦ 

Subscript 

c 

CG 

u 


Vorticity strength per unit length or ratio of specific heats 

Horseshoe vortex strength in Trefftz plane 

Control surface deflection 

Kroneker delta 0 i^j 
1 i=j 

Jet deflection angle relative to trailing edge 
Total jet deflection angle 

Body slope 

Arc of jet segment 

Dihedral angle of quadrilateral panel or boundary layer 
momentum thickness 

Absolute viscosity 

Kinematic viscosity, y/p 

Density 

Source density 

Perturbation velocity potential 
Total velocity potential 
See figure 3 

camber 

center of gravity 
upper surface 
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list of symbols (CONTINUED) 


I 

LE 

r 

t 

TRAN 

00 

Superscript 

I 


lower surface 
leading edge 
recovery 
thiclciess 
transition point 
freestreara condition 


quantity based on effective origin 
Eckert reference ten?>erature condition 
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THEORY 


The arbitrary configurations which may be treated by the analysis 
are simulated by a distribution of source and Vortex singularities. 

Each of these singularities satisfies the linearized small perturbation 
potential equation of motion 

A 4> * d> 

The singularity strengths are obtained by satisfying the condition 
that the flow is tangent to the local surface : 


All of the resulting velocities and pressures throughout the flow may be 
obtained when the singularity strengths are known. A configuration 
is cdn5)osed of bodies, interference shells and a.erodynamic surf aces (wings , 

canards, tails etc,)^ The following types of singularites are used to 
represent each. 


wing and vertical tail 



The first step in the solution procedure consists of obtaining the 
strengths of the singularities simulating the fuselage and nacelles, from an 
isolated body solution. The present analysis uses slender body theory to 
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predict the surface and near field properties. The solution is con 5 )osed 
of a compressible axisymmetric component for a body of revolution of the 
same crossectional area and an incon^ressible crossflow con^jonent, cj) , 
satisfying the local three dimensional boundary conditions in the Cy,z) 
plane. The crossflow is a solution of Laplace *s equation 

+ <P S! O 

A two-dimensional surface source distribution fomulation is used to 
obtain this solution. When the body singularity strengths are determined, 
the perturbation velocities which they induce on the aerocfynamic surfaces, 
or other regions of the field, are evaluated. 

pie assumptions of thin airfoil theory allow the effects of thickness 

lift on aerodynamic surfaces to be considered independently. Therefore 
the effects of the aerodynamic surfaces can be simulated by source and 
vortex singularities accounting for the effects of thickness and lift 
respectively. The source and vortex distributions used in this program 
are in the form of quadrilateral panels having a constant source or vortex 
strength. The vortex panels have- a system of trailing vorticies extending 
undeflected to downstream infinity. The use of a chordwise linearly 
varying source panel is provided as an option to eliminate singularities 
associated with sonic panel edges at supersonic Mach numbers . 

The panels are planar, that is they have no incidence to the free stream 
(although dihedral may be included) , since thin airfoil theory allows the 
transfer of the singularities and boundary conditions to the plane of the 
me^ chord. These boundary conditions are satisfied at a single control 
point on eadi panel. For thickness, the control point is located at the panel 
centroid while the effects of twist, camber, and angle of attack are . 
satisfied at the spanwise centroid of each vortex panel and at 87.5 
percent of its chord. 

A cylindrical, non-circular, interference shell, composed entirely of 
vortex panels, is used to account for the interference effects of the 
aerodynamic surfaces on the fuselage and nacelles. The boundaiy conditions 
on an interference shell are such that the velocity normal to the shell 
uiduced by all singularities , except those ^pf the body which it surrounds, 
is zero. The boundary conditions are satisfied at the usual control points 
for vortex panels. 

The following sections define the details of the solution procedure. 
Included are discussions of the isolated body analysis, panel geometry, 
boundary conditions, and influence equations, the jet flap solution, and 
evaluation of aerocfynamic characteristics including drag. References are 
cited for the reader interested in further pursuing a particular point. 
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Body Solution 


According to slender body theory^? 2 -thg flow disturbance near a 
regular three dimensional body may be represented by a 
perturbation potential of the fom 

+ 5Cx) 

is a solution of the 2-D Laplace equation in the y, z cross flow 
plane satisfying the following boundary conditions 


V<p = jv + = 0 

3 $ 

‘ ° on C(x) (2) 

C(x) and n, are defined in figure 1 . A general solution for <f may be 
written as the real part of a complex potential function W(Z) with Z * y + iz. 

<P « W « r^c{ A.cx) JU.Z A_cx)z’"^l 

A useful alternative representation of and W is obtainable with the aid 
of Green’s theorem? 

<}> ' - 2 fuj JU(Z-l) (3) 

CCx> 

where <r(s) is a ’’source” density for values of 1 * y^ + iz^-, (yc*Zc) being 
coordinates of a point on the contour c(x). 


The function g(x) is obtained by matching <p of eqioation (1) which is 
valid in the neighborhood of the body with an appropriate ’’outer” solution. 
gCx) is then found to depend explicitly on the Mach number M and 
longitudinal variation of cross sectional areas S(x) 






TT { J -h f 5"ct) JU^ct-x>at 

X 

- J 5'co) JU.X - ~ 5C0 j 

X 

a cx) ) - J 3'*ct> jlv^.Cx-t) cAt I 


M < I 


C4} 


n > I 


where 


/b . Vl I 
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The body axis perturbation velocities are obtained by differentiation 
of equation (1) 


V ^ i> 

lo , <p^ 


At supersonic speeds , -zone of influence considerations require that u » v =* 
w = 0 for x-ydt* < o. 

Solution of the preceding equations is based on an extension of the 
method of reference k . 


Cross Flow Component 

The reduction of computations to a numerical procedure utilizes the 
integral representation of <p given in equation (3) by. discretization of 
the cross sectional boundary into a large number of short linear segments 
(figure 2) over each of which the source density <r is as.sumed constant at 
a value determined by boundary conditions. 

Computation of <T(i,n) over the segment i, i+1 proceeds by applying 
the boundary condition equation ( 2 ) at each segment of Cn. If vf =q*jv+kw 
represents the velocity vector, the corresponding con 5 )lex velocity in the 
cross flow plane is obtained by differentiation of W in equation (3) with 
respect to Z: 

v-lu, . - (S3 


The contribution by the sources located on segment i, i+1 to the velocity 
,n first evaluated. Noting that i, i+1 makes an angle 0(i,n) 
with respect to the horizontal axis, we Imve 




o\.A <2. 


and the contribution to the integral in equation ( 5 ) may be written: 






1. 


Z. - i 

4,«n, 
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^n+l ^n-l 


Figure 2, Crossection Boundary Segmenting Scheme 
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After integration of the last tern and summation over all contributing 
segments, the result may be written 








in which, referring to figure 3 , the quantities R(i,j,n) and g(i,j,n) 
are defined by the relationships 

e = Z . - i. 


To insure uniqueness of the coiT5)lex velocity, care must be exercised 
in assigning values to the angles 'f(i,j,n) and rU,j,n). Referring to 
figure 3 , these are measured counter-clockwise from the positive y 

axis so that when facing from ^ to Pi+i,n , a point ^ 3^^ to the 
left of i,i+l shall define an angle 'f'(i,j|n) =* ^(i,n). *As Pj 
traverses a path around P^^n to a point just to the right of i,i+l, V'(i,j,n) 
increases from 5 Ci,n) to d(i,n) + 2 ir . The same holds true for f'(i,j,n) 

as Pj,n traverses a path around Pi+i^n* consequence of these definitions 
SCi,j,n) becomes -tt when approaching i,i+l from the right and 77 when 
approaching from the left. This discontinuity reflects that exhibited by 
the stream function upon traversing any closed path which encloses a 
distribution of finite sources. 

From the boundary condition equation ( 2 ) , we have 





After substitution of v and w from equation ( 6 ) , this last expression 
becomes 


where 




4. -vs. 




( 7 ) 




+ [ 0(4,'>v) - ^ 
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Figure 3. Details of Variables Pertaining to Segment i,i+l 
of Boundary 
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The surface normal perturbation velocity - may be written in terms 

of the body slope CaV/ax)j^^ , the angles of attack , and sideslip a 
and the angular velocities p,q,r as * 



PjjJ 6(4,^:) 




Satisfying equation 7 at each of the points P- ^ on a given contour 
boundary yields a set of equations for cr (i ,n) . ’ 


Axisymmetric Component 

Differentiation of g(x) must be carried out with due concern for the 
nature of the improper integrals appearing in equation (4) . The result is 


where 


3'<-0 = I s"c>Q ^ - J^OQ 

- s'<0 )* 5"coJU(I-*o} f1 

^ ■* I ■^^”*•0 ~ - ■ s ”< o 4 JLvx^I m 


M-i 


] JUCx^-t)5'7t)dt = 

y [ s'- s^i X •) 




I 

0 

2 [ Ct 

-m. : 0 

* T ^ ^ ) 
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To compute the second derivatives of the equivalent body cross sectional 
area required for g^(x),the first derivatives at are found by finite 
differences between and %+x« Second derivatives at x’k 

(xm+l + are then found by finite differences between S’ at xji and 

, Finally S”(%) is determined by linear interpolation of S”(x^) 
between and 

Perturbation Velocities 

The axial velocity u depends on (af/ax) and the axisymmetric solution 
g’^Cx). (af/ax) is obtained by differentiation of the integral in equation 
C 3 ) to first obtain an exact expression which is then approximated by 
evaluating the result over the segmented boundary. 


The derivation of a<P/ax must take into account the fact that the 
path of integration in equation (3) is a function of x. Referring to 
figure 1 increments of a dependent variable taken along CCx) are denoted 
by d( ) and incronents taken normal to C are. denoted by ). Differentia 
tion of equation (3) then yields 


■|r * ■ ^ i If - j 

From figure 1 

= S-V^e = St' 4^, C93 


where h(l) is the radius of curvature of C(x) at ^ . In addition, we have 
from figure 1 


Sx " Sx (10) 
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To evaluate -^-21 we note, 
Sx 


ISL 

8x 


8x-»o 


t ) - <T( i 

Sx 


(11) 


Introducing equations ( 9 ) , ( 10 ) , and ( 11 ) into equation ( 8 ), 



~ I 7 ] 



Again, assuming that quantities in the brackets of the integrands are 
constant over i,i+I, 






vdiere 


- ^ S ( ^,4 




■| R ( J 4 i'W) • U c t, ) >^/vs. R( t+!, 4 ,-»v) 

RC*-,4,*n.) •Us(tj-n.) JL^ R.C^>4,'>^) 

- RU,4,nx) SCi,4,-vv> 
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The radius of curvature h(i,n) and the derivatives S«*/Sx , Sv/gx are 
approximated at the raid points of the segments i,i+l as follows 

a) S’T/ Sx - the derivative at the niid-point x'n of the interval 
^n,^n+l is set equal to the divided difference between <r(i,n) and ^s'Cijn+l). 
Linear interpolation between these derivatives then yields 6 <r/ 5 x 

at Xn* 

b) Sv/Sx - referring to figure k , the displacement $7 is 
determined by linear inteipolation between S 5 ^ ^ and 5 5 i+i n- 

^J/O^n+1 - Xn) then represents Si/y^xat x’n* ’Linear interpolation between 
the stations x»n then yields Sv/sx at x^ . 

c) 1 /h: - 0 at is determined by interpolation between values of 
y(i,n) at P’i,n- The curvature 1/h at P»i „ is then set equal to the 

divided difference between & at P^+^^ ^ and & at Pi^n . 

The lateral and vertical perturbation velocities, v and co are 
obtained from * 


V - Z lO 


- a i ^ 
J z-l 


Integration over the boundary with constant segment source density yields: 




L J 


Thus 



z 










lO ■» 



ft 


RC**I,4 




^ ( t ,4 ,-w) Cfr<» 
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Panel Singularities 


The source panels and vortex panels are composed of quadrilaterals 
with two edges parallel to the free stream. The coordinates of these panel 
comers are specified with respect to an (x,y,z) system having its x axis 
in the free stream direction and its z axis in the lift direction. However, 
the panel influence equations are written in terms of a coordinate system 
having a z axis normal to the panel and an x axis along one of the two 
parallel edges. A coordinate transformation is necessary to obtain the 
coordinates in the panel reference systam. If the plane of the panel is 
inclined at an angle 6p with respect to the y, z plane, a transformation 
into the panel coordinate system (x^,yp,Zp) is accomplished as follows: 


^ ^ Co^ * 2 




control 

point 

panel 





A transfoimation of the (t^,Vp,v^) velocities into the coordinate system of 
the panel on vMch the control point is located (uc,Vc,Wc) results in the 
axial, binomial and normal velocities induced on the panel. 

For the image of the influencing panel, the signs of y, and v^ are 
changed while using the same calculation procedure. 
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Panel Singularity Strengths 


The source singularity strengths may be foimd directly by equating 
each source panel strength to the slope of the thickness distribution at 
its control point. For panel i 



where refers to the shape of the thickness distribution. The influence 
eqiiations for the source panels can then be used to obtain the velocities 
induced by the source panels anywhere in the flow. 

The determination of the vortex panel singularity strengths are the 
final step in the solution procedure. They are obtained by solving a set 
of simultaneous equations utilizing the vortex panel influence equations to 
relate the singularity strengths to the boundary conditions at the control 
points of the vortex panels. The boundary conditions permit the condition 
of tangential flow to be satisfied. 


Each vortex panel j having singularity strength Cp. induces a set 
of velocities ( , A^- , A?. ) on panel i. Therefore a set of influence 

equations can be written: * 


-c, - 2 A,; c,. 

A 4 4 

+ u. 

o. 

"s • f Ar, 


4 '* 

^ 0. 

i> 

where (uoi* Voi, wq^^) refer to the velocities induced by all other body 
and source singularities, and written in the coordinate system of the panel 


containing the control point. Since the resultant velocity along the normal 
at a panel control point must be zero. 
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and the following system of equations results 






A . . c. p 

<•4 4 



This set of linear equations can be solved for the Gp. and, since it 
assumes symmetrical panel loading, can be used to determine the longitudinal 
characteristics. A similar set of equations exist for the calculation of 
the lateral/directional characteristics. This set assumes an antisymmetrical 
panel loading and has a correspondingly different set of influence 
coefficients Aij . 


Boundary Conditions 


Several types of basic and unit boundary conditions are considered 
and can be classified as either symmetric or antisymmetric. Linearized 
theory allows the superposition of these basic unit solutions. The p, q 
and r rotary derivative boundary conditions are the result of placing the 
configuration at®c=0, y6=0ina flow field rotating at one radian 
per second. 

Symmetric: 

1) basic 

S *a 


Cdls^cLx') 


« surface slope due to twist and camber 



» noimalwash induced by slender 
body thickness and camber 


• noimalwash induced by source panels 



2) IMit alpha 


3) Ifeit q rotation 

4) Iftiit flap 

Antisymmetric : 

1) litiit beta 

2) IMit p rotation 

3) lAiit r rotation 

4) Ifeit flap 


- CO 

ISO ^ ® 

» normalwash induced by slender body 
at unit alpha 

- -=r(K')C ) c-oA© - (O 

c 

» normalwash induced by slender body 
undergoing unit q rotation 

TT 

- <r 

Igo 

^ » 1. for flap panel 
<T“ =0. for others 


■nr « ^ 

~ *■ Q — A 

i8o ^ 

» normalwash induced by slender bocfy- 
at unit sideslip 

- C 

» noimalwash induced by slender body 
undergoing unit p rotation 


^ Cx-Xe^) u> 

- Normalw^h induced by siajder body 
undergoing unit r rotation 

TT ^ 

— - cr 

l6o 

^ * 1. for flap panel 
^ ■ 0. for others 


23 



Constant Source and Constant Vorticity Panel Influence Equations 

A constant pressure or constant source panel with a quadrilateral shape 
can be constructed by adding or subtracting four semi -infinite triangular 
shaped panels . These semi -infinite triangles, each determined by a comer 
of the quadrilateral, can be assumed to induce a velocity perturbation every- 
where in the flow. However, each comer represents only an integration limit, 
and all four comers must be included to make any sense. 



If it is kept in mind that foxir comers must be included, one of these 
triangles having sides determined by y * 0 and x-Ty = 0, induces the 
following perturbation velocities: 
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I /3 > o 
L 2 /i < o 

S I - M*- 

K s x' + + z*) 

Constant source panel 


w-Cx,^,E,T) 






^<r » /rVya* R 4. 

2ir 2 


^ir 


it<r 


^ oVt^*- -(Tx+/a'3^ 


^1 
'3^ J 


-I 


z R. 


^ *3 -fCaV*') 


Constant vorticity panel 




Vci 


Bit 


{ 




-I ZR 






■f} 


Vt>^,3,4^T) = - 




^[tX^-' £5. 

8'^ I 


•^C|i 

87T 


7-r + 4ri-^--*‘>fT:^"-| --f^li 

>€*) L 2 

I ^ ^ R-x e R - CTx»/i*3) 


25 



These perturbation velocities hold for both supersonic and subsonic free 
stream velocities. In supersonic flow only the real and downstream contributions 
are considered. 


To establish that these are the correct perturbation velocities the 
following criteria must be met: 


1. Laplace's equation must be satisfied 




or the equivalent. 


2€ 


u. < V 


u., 


lO. 


V, » tO 


u., . y ♦ 


2. The correct discontinuity or jump in the perturbation velocity must 
occur at the surface of the quadrilateral panel area. For the source 
panel the jump occurs in the normal or w velocity and on the vortex 
panel there must be a jump of constant magnitude in the u perturbation 
velocity over the panel area. The perturbation velocities should be 
continuous elsewhere, except on the trailing vortex sheet of the 
vortex panel. 


3. The perturbation velocities must go to zero as upstream infinity is 
approached . 


4, For the vortex panel the trailing vorticity must extend straight back 
to downstream infinity. This means that any discontinuity in the v 
velocity must be zero outside the spanwise boundaries of the panel and 
must be zero upstream of the panel. 
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The first criteria can be established by using the derivatives given in 
TABLE I. 


The second criteria can be establisiied by noting that all terms except 




-I 


— and 


are cont^uous at Z G. Consider these terns keeping in mind that the 
contributions from all four comers must included. 


If we let 


C.JC-X,) - T( , C x-x^) - T ) 

' Cx-X.) + /ft' t C + l'] 

and use 

^^■'a . Wb 3 x >^' 

I - Ae 


then the contributions from both comers on the leading edge can be combined 
as follows. 




- ra*- 




X -Ta' 





t - te‘] R, - ) -Ta‘3 

— : . 

I - T 2'lt 5 c V'1,> - ra' j 2 *r,R 
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TABLE I 


DERIVATIVES 

R.* = it'+ 



-L J. Jt- JLL2L 
ax i o R-x 

_a ! i_ R»(Tx*jS*a) 

3x VtV^‘ 2 ^v'r^n-Crxt/a’a) 

a + -1 HR. 

3X 


R 

H(Tx-»/a*a^ _L 


R.- 

R = 





d 


R + X 



\ 


z 

R-x 



) R 

d 

1 

1 0 *ASya^ R 

4. CTxi->s*a> 

[x-U-Ta>+^V] , 

»3 


2 ^ R 

- CTx-^/a':^) 

cx-Ta)* 

* a* R 

d 


HR 


X z 1 

^ TECTx+/»’a) 1 




;i^**a* R 

tX-Ta>'+ (^’* t‘)2* r 

R 

^ R = 
9i 

/3*H 




d 

— JL«. 

RtX 


xa 

1 

d£ 

2 "o 

R-K 



R 

d 

1 

, 7?^ R + 

1 

X 

HCTx-^A\) 1 

d£ 

V^V 

2 ^Vt V R - 

<Tx*;a’a) 

Cx-Ta>' 

♦ Cyft'^T*) z*- R 

a 

_ 

HR 


1 

<J^"T:iKTx4/a'3) ^ 
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If we define 




0® -4. -I a. 

Aol^ — 



- 7 T i c i TT 


A-i (E) 


■” 4*C4) 

z -^o* a -*o’ 


then 


SIZ. 

Z-^0 


{ ju- 


Z R, 






2R, 


scvv*^* 


■1 




o 


(-TTa^vZ 




**• - 

w * '*• ] 

o 

r 


[ ICv^.)-T£» 


- Zrr 



o 



Therefore when a similar procedure is carried out for the trailing edge of 
a source panel we obtain the following jump in the w perturbation velocity. 
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For the vortex panel (subsonic) we have an additional term 
both additional terms from the leading edge corners j 


Cons idering 
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Therefore combining the terms 
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The contribution from each panel comer is : 


A tv. 
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8 tt 
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Therefore summing all four panel comers 
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To verify the third criteria we must show that all of the functions 
approach zero when all four comers are considered as x — « 









Therefore considering both comers on the leading edge of the panel 


R *» ixl 

X - o- 




X - •• 



R, » Cx-x.^ 
R, 4 (x-x, ) 





Cx-x,^ 

Cx-x^) 





R, i> Cx-x,^ 


o 


Of,,. 
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and therefore this limit is also zero when both comers of the leading or 
trailing edges are considered. Since all terms are accounted for, the 
perturbation velocities are zero far upstream. 
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Since 




« >/t\^ r. + (tx + /s':i) 

^ A *■"/»•' - CTx+/i';i) 


^ /»* t C A -Ta)\ (T Vya' ) a' 3 

A^ [ U -Ta Ct\A^) 2^ 1 

r< - 



there is an apparent singularity along the line 

» o ^ a » o 


However this singularity may be 
both corners of the leading or 
these edges the values of 


removed by combining the contributions from 
trailing edges of the panel. Along either of 


are the same for each of the panel corners. 



It can be seen from the above diagram that tr» +/i*a) the sen,e 

auaLi?^r outside the spanwiae boundaries of the ^ 

quadrilateral. Therefore outside the spanvise boundaries the term 

can be eanceled by combining both comers, and the resulting term 



R. t CtCx^x.) ■> 
t 


will not be singular if the correct + or - sign is chosen. Within the spanwise 
boundary an actual singularity occurs on the panel edge. 


33 



The term J log 
can be written 


also has a possible singularity. This term 


ilog 


R*X 

R-x 



(Ri-x)* 


For the source panel the singularity may be removed for points along + 
which are outside of the panel boundaires. 


If ( x-x,) and ( x - ) have the same sign the combination of the two terms 

gives 


^ R, -Cx-x.) “ I 


Cx-xQ 

-Cx-xO 





R, t Cx«x,^ 
Rj 1 Cx-x^) 


where the correct sign is chosen to remove the singularity. On the panel 
edge the singularity is real and cannot be removed. 
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removable singularity 
real singularity 


For a vortex panel the terms (subsonic) 


3 C R.»x) 


Both have real singularities for x>o (downstream) and removable singularities 
for x<0 (upstream). The real singularities occur on the panel edges and on 
the edge of the trailing vortex sheet. 
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Supersonic Velocities - Special Considerations 


The velocity perturbation influence equations for supersonic flows are 
treated by taking o nly the real parts of the expressions. This means that 
^ = V X is set equal to zero for points which lie outside the 

downstream njach cone from any given comer. Therefore R and J log a- 

zero for points which lie outside the downstream Mach cone. For Ct*- >o 
there ere no problems using this method. ^ 


li* O* “ ® the real part of 
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therefore combining two corners 
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Tf a - 0 and either R, or R is zero and we allow the other to approach zero 
the value of F2 becomes 


F2 - 




Therefore if R^ and R^ 
from the 3eading edge 


are 

(see 


zero but we are inside the envelope of mach cones 
figure 5 ) the value of P2 is set equal to 


F2 


TT 

7^ 


if 


[tu-x.) - 

-• R, < o < o 

> C>a‘-T‘) a* 


< o 


35 



sdoxsAUH 3003 sSpg 3irtp33X oxuosasdis S sanSxd 




< lOtL-tia-fK-o] . V— 5f -X 
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Linearly Varying Source Panel Influence Equations 


In st^Dersonic flow constant source panels having a sonic edge have a 
real singularity along an extension of this edge. The singularity occurs 
because : 

-» o 




Control points which are near the extension of this edge will have large 
u and V velocities induced upon them. The singularity can be eliminated by 
using panels vMch have a source distribution \diich varies linearly in the 
chordwise direction, llie resulting continuous source distribution eliminates 
the singularities. The linearly varying source panel influence equations can 
be found by integrating the constant source panel influence equations with 
respect to x. 
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These velocity con^jonents satisfy the same criteria as the velocity con?)onents 
for the constant source panels except that the source strength is proportional to 
x-Ty. The source panel finite elements are constructed with the following 
properties . 

1. All panel leading and trailing edges are at constant (^) , side edges 
are at constant y. 

2. Each source finite element is con^josed of a pair of chordwise adjacent 
panels . 

3. The source strength varies linearly with chord measured from the 
leading edge of a panel pair, i.e. the maximum value of the source 
strength is proportional to the local chord and attains this maximum 
on the panel edge joining the panel pair. 



The perturbation velocities induced by this panel pair are conmosed of 
contributions from six comers. 




V3,,*,T*) - w-,.U-X3, 
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If there are N panels in the chordwise direction there will be N-1 
singularities or unknown source strengths associated with them. The 
linear variation in the source distribution means the value of dz/dx 
must be zero at the leading and trailing edges of each span station. 
This may be an undesired restriction and therefore the use of linearly 
varying source panels is optional. 
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Numerical Solution 


Householder’s method for solving simultaneous equations is used in the 
solution of the aerodynamic influence equations . The influence matrix is 
triangularized by means of orthogonal transformation matrices, which 
preserve the conditioning of the matrix. This along with a reduction in the 
number of required con^juter operations greatly improves the numerical 
accuracy and stability of the solution over that of the standard Gaussian 
reduction method. 

A complete derivation of the method is given in reference (6) . The 
method has been altered' from the original to allow the operation on a 
single row of the matrix at a time. This reduces the required core allocation 
necessary to triangularize the matrix. 

If [A] is the square influence matrix, the upper triangle is given by 

fR] “ CW] [AI 

where [W] is the combined orthogonal transfoimation matrix used by House- 
holder to triangularize [A]. 

In the Householder method [W] is equal to the product of N individual 
orthogonal transformation matrices, where N equals the number of unknowns. 

Each transfoimation results in reducing all elements below the diagonal to 
zero for one column. The columns are reduced from left to right. 

The individual transfoimation matrices [W]m are defined by 

(fl] - 2{u„,){%r) 

where [I] is a imit diagon^ matrix and iHml is a column matrix defined by 
the unit vector % (a^ Jijjj.xhe vector % is defined by the mth 

column of [A] where the elements on rows less than m are replaced by zeros. 

The unit vector ^ is defined by a column matrix |Vjjjl with all zeros 
except for the mth row, which is equal to one. The constants and Um 
are defined as 

«m “ I l%l I 
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It can be shown that|anj)is reduced to ]{ani}( {vm} if {am) is premulti- 
plied by C [I] - 2 Also, that the first m-1 rows of 


[Wm-l][Wm-2] ••• WEA] 

remain unchanged by the mth transformation. The result after m transforma- 
tions is then zeros below the diagonal for the first m columns and 
|{a}il , [ {a> 2 [ . , [{a}m-2| , [{a}m-l( , |{a}m| on the diagonal. The 

elements above the diagonal have been defined by the m preceding 
transformations and will remain unchanged for the N^m remaining 
transformations , 

The relation 


(UI - 2|unili%F) 1%! “ ii%inv 


( 12 ) 


where - or^VjO/Mm or |%1 = ( 1 %) - an,jVj,|)/ 

“in ” 

and 

*‘m * ' %•%) or - '^2o-^(a^ - |ajj)T|v^|) 


remain to be proved. It is helpful in the derivation of equation (12) if 
the vector identity 

|^|V+ 20^.1;^ 1^ = ^ 


is observed from the following vector diagram. 
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Then from equation (13) 


l\l^m * 2 ^Cum • ^) = \ 

f^l^m * 2 ^ = a;„ 

Therefore 

(143 

where and i%i%n are dyadics . The unit vector is in the direction 
of aJn* 

Equation (14) can then be written in . matrix notation as follows 


([I] - 2 |%|(Umr)(aB,l = llanlHv^) 


which is equal to equation (12) . In matrix or tensor notation it becomes 
evident that the dimensions of {a^}, {Vjj^}, and {i%} are not limited to 
three . 


“m “ IKII 

and 


Then if equation (12) is premultiplied by {am}^ 


lamFlaj,) ' 2130^1% 
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And substituting and into equation (15) 


»'m= '^2am(aa, - 


In vector notation equation (16) is seen to be equal to 
Also, if equation (16) is substituted back into equation (12) 


1^1 - - |aj,|T|vn,l) luj = 


Therefore 

, , - “mNI 

t“ml = . 


or in vector notation 



Jet Flap 


A completely linearized approach was used based on the assumptions of 
thin airfoil theory. The flow was assumed to be inviscid and irrotational 
and all entrainment effects were neglected. The jet was represented by an 
infinitesimally thin sheet having zero mass flow but finite momentum per 
unit of span. This sheet was assumed to extend from the trailing edge of 
the surface back to infinity. (In practice one or two chord lengths is 
sufficient) . The effects of transverse momentum and the deflection of the 
jet sheet were neglected. 

Since both the planform and jet can maintain a pressure discontinuity, 
they are both represented by a system of quadrilateral panels having 
continuous distributions of vorticity. The strengths of these constant 
pressure vortex panels are determined by solving a set of linear 
simultaneous equations which satisfy the downwash boundary conditions at 
a set of control points on the planform and jet. 

The boundary condition on the planform is the previously described 
flow tangency condition. The pressure difference across the jet causes a 
change in the direction of the jet momentum. The equation relating these 
quantities forms the boundary condition on the jet and can be derived by 
considering a jet segment of unit depth. 



The mass rate of flow through the jet is m and the velocity is v. If we 
assume a pressure difference of AP across the jet, then from the momentum 
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theorem applied to the differential 


m 


or 


VA54 ^ 


element, we write 
^PR Ajrf 


= 


■mV 

R 


where R is the radiiis of curvature of the jet. 

The reaction of the jet on the flow external to the jet is 

F = aP RAji 


A vortex of strength per unit laigth along the jet of y would produce a 
reaction of 




Hence, equating these two forces, we calculate the action of the jet on 

the flow external to the jet by replacing the jet with a running vortex 
strength ^ given by 


/>UJ< 

For a nearly horizontal jet with a large radius of curvature 


1 w 

R ^ d X2* = ax 


where w is local downwash velocity (nondimensionalized with respect to lO 
Then 


± C 
2 
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(Y) C (Y) ( X,Y) - ‘^Pr*ET (X,Yj = O 

which is the boundary condition written for a three dimensional jet flap. 

To apply the jet flap botmdary condition to control point i, the above 
equation is integrated between adajacent control points in the streainwise 
direction. 

9,(yJc(Yjj l^(x,Y)ax-| '^Pnet (x,y) ax - o ( 17 ) 


The control point is located at 87.5 percent of each panel chord. 

To siiiplify the second integral in equation (17) the assun^jtion is made 
that the control point is exactly at the panel trailing edge. The effects 
of such an assuti5)tion have been shown to be negligible. Equation (17) 
evaluated from the leading to the trailing edge of panel i yields the 
following relation: 









o 


(18) 


The downwash at each control point is written in terms of the N net 
presstires on the quadrilateral panels: 

N 


fuj 


1=^ AXj 
i = \ 


AJlT 


n1 

'i-' 

Equation (18) is then written 


N 


where is the Kroneker delta. 


47 



For a flap panel adjacent to the jet exit, equation (18) must include 
any jet deflection angle relative to the surface trailing edge. 

-(rL-i ^-9^)] - AX^cO 

or 

C^C ^Xij ^PNETi = 

where is the jet deflection angle. 

Then 

E C.C (K j -jl - AXl] c 

i -1 


The complete set of linear simultaneous equations for both the surface 
and the jet flap is then written 


Vi 


^PlsiET^ = 


1, M 


where 


FI 4. 


C«C [Aij-Ai.-! j] AXL for i on the jet 


and C-^ = 


A 4 

for 

_c^c 


C- <^i 

for 

9 -cr, 

for 

L 0 

for 


(19) 


Both syinnetric and antisymmetric jet deflections are considered. Thus 
after calculating the influence matrices and boundary conditions in the 
usual manner, the appropriate rows are modified and combined to produce a linear 
symmetric or antisymmetric system as described by equation ( 1 ^ . Because 
of the rotational quality of the flow fields, the p, q and r rotary 
derivative calculations are generally not valid for jet flap configurations. 
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AERODYNAMIC CHARACTERISTICS 


Longitudinal and lateral -directional forces and moments due to thick- 
ness, twist and camber, pitch, sideslip, and the dimensionless rotary- 
velocities are obtained from surface pressure integrations of the 

various configuration components . 


Bodies 

The pressure coefficient, to an approximation consistant with slender 
body theory, is 
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( 20 ) 


The forces and moments are obtained from the surface integrations 
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are 


In terms o£ these expressions, the commonly used aerodynamic coefficients 
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where L is the body length and c ^ b and S^g£ are configuration reference 
chord, span and area respectively. 


50 



CrosscoL5>lii^g between the pitch, sideslip, rotary motions through 
> the product and quadratic tenns in equation (20) is neglected. 


Planar Components 


Surface pressure distributions are calculated for planar components 
using the first -order linearized form of the pressure coefficient. 





* 


The^V* signs refer to the upper and lower surfaces respectively. The term 
consists of the velocities induced by the isolated bodies and 
other vortex and source panels. These velocities are obtained by multiply- 
ing the u" influence matrices by the appropriate panel strengths.' fhe 
term accounts for the ■§• perturbation velocity induced by the 
local distribution of vorticity and changes sign from upper to lower 
surface. The total and values are the result of taking linear 

combinations of all the basic and unit solutions. 

The net pressures for each of the basic and unit solutions are integrated 
numerically to give the section forces and moments, component forces and 
moments and configuration forces and moments. 

Since the vortex panels have a constant pressure distribution, a block 
integration scheme is en^Dloyed. With the exception of drag, these basic 
and unit force and moment coefficients are combined in a linear manner to 
produce the aerotfynamic characteristics for any desired flight condition. 

Since drag varies in a parabolic manner, it must be considered on a point 
by point basis as defined in a later section. 

The longitudinal normal force distribution on the bodies in calculated 
for each solution. The load distribution on the interference shell portion 
of the body is given by integrating over all vortex panels at a given 
longitudinal station. 

normal force 


T * AX L X! 
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where n is the number of panels around the shell, L is the length of the 
body, AX is the length of the interference shell segment. At is the panel 
area and Ci = 2 for a centerline body or = 1 for an off centerline body, 

piis carryover load distribution is added to the previously calculated 
isolated body longitudinal load distribution. 

The section characteristics of planar con5)onents are determined by a 
chordwise sunmation of panel data at each span station and are given by 
the following equations: 

local lift coefficient 

Ji^ 

C = — ^ A' 

CAS / 

1=1 

normal force 

N 


lift force 


'P 7^ = V e; 

C>\v6 As y 


AM 


center of pressure 


C AsC 




AVCi 


N 

E 


A- I 


^net: Ai ) 


•n 


“AVG 


where N is the number of chordwise panels, and As is the width of the span 
station and is given by 

AS 
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The section characteristics due to the reaction o£ a jet flap are 
calculated by taking the appropriate component of the reaction force. 


reaction normal force 


■n 


Cavg 


= 9- 


Cavg 


Sin 


c 






where is the total deflection angle of the jet. 


reaction lift force 




- C. 


Cavg 


COS ©- 


Component forces and moments (excluding jet reaction loads) are given 
by the following equations: 


lift 
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L 



Ax COS 


side force 
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rolling moment 


N 

1=1 

pitching moment 

. N 

^ g 5r£f-^ ' ^Cg) 

i = l 

yawing moment 

N 

<^A = ^ ‘^PneT; Ki SiN (x • - Xcg) 

1= I 


where N is the total number of vortex panels on the component, and F^ and 
F 2 are defined as follows: 

symmetric loading F^ = 1 asymmetric geometry 

= 2 symmetric geometry 
F2 ” 1 asymmetric geometry 
= 0 symmetric geometry 
antisymmetric loading F-j^ “ 1 asymmetric geometry 

= 0 symmetric geometry 
F 2 = 1 asymmetric geometry 
- 2 symmetric geometry 
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The X coordinate of the center of pressure is given by 


C 


^C.R - 


XG 


For interference shell car^jonents , the total forces and moments of the 
corresponding isolated body are added to those of the shell. 


Jet reaction forces and momfflits are obtained from a spanwise 
of the jet flap section characteristics: 


summation 


lift 
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4.= I 
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yawing moment 



where N is the number of spanwise jet flap stations and Fx and F2 are 
defined as before. 

The forces and moments for the complete configuration are obtained 
by summing those of the individual components. 
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DRAG ANALYSIS 


Preliminaiy estimation o£ configuration aerodynamic efficiency requires 
the calculation of aerodynamic drag. The program separates the conroutation 
into two parts which are assumed to be independent of each other: viscous 

skin friction drag and pressure drag. The .latter calculation results in a 
set of drag polars representing the attainment of zero or one hundred 
percent of the leading edge suction available. The conplete estimation* of 
the pressure drag requires an assessment of the percent of leading edge 
suction which can be attained. Although the loss of "suction” is due to 
viscous effects, it is assumed the suction level will be estimated semi- 
enpirically and therefore is not part of the viscous calculation. Trim drag 
polars may be calculated to assess the tradeoff between longitudinal 
stability and performance. 


The calculation procedure assumes the viscous skin friction drag and 
supersonic pressure drag d*je to thickness are independent of lift. The use 
of a linearized analysis results in a quadratic drag polar for the 

drag due to lift. Therefore. , the trimmed or untrimrued drag polar for 
zero suction or one hundred percent suction is of the form 



Co + C* + Cb > 

TMickm«4s Mt«» 


KCCu-Ct )* 


where K depends on the type of polar. The specific techniques used to 
acconqjlish the drag analysis are discussed below. 


Skin Friction 


Several well established sonioipirical techniques for the evaluation 
of adiabatic laminar and turbulent flat plate skin friction at incompressible 
and compressible speeds are used to estimate the viscous drag of advanced 
aircraft using a con^onent buildi^ approach. A specified transition point 
calculation option is provided for in conjunction with a matching of the 
momentum thickness to link the two boundary layer states. For the 
turbulent condition, the increase in drag due to distributed surface 
roughness is treated using unformly distributed sand grain results. 

Conqjonent thickness effects are approximated using experimental data 
correlations for two-dimensional airfoil sections and bodies of revolution. 


Considerations such as separation, component interference, and discrete 
protuberances (e.g. antennas, drains, aft facing steps, etc.) must be 
accounted for separately if present. 
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In the following, a discussion is presented for a single coiiponent 
evaluation in order to simplify writing of the equations and eliminate 
subscripting. The total result is obtained by a surface area 
weighted summation of the various conponent analyses as described on 
page 63, 

Laminar/Transition 


A specified transition option is provided in the program. The principal 
function of the calculation is to provide the conditions required to 
initialize the turbulent solution. In particular , the transition point 
length and momentum thickness Reynolds numbers are required. 




R — ^ L 
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where 
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T+ 1^8.6 


lb sec/ft^ 


This solution is based on the laminar Blasius result C8, chapter VII) 
in conjunction with Eckert *s compressibility transformation^. This 
option pemits an assessment of the reduction in skin friction drag if 
laminar flow can be maintained for the specified extent. It does not 
establish the liklihood that such a condition will be realized in practice 
or to what extent. 
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Turbulent 


Smooth and distributed rough surface options have been provided in the 
analysis. In either case, the solution is initialized by matching the 
momentum thickness at the transition point produced by the laminar/transition 
solution. That is, an effective origin (commonly referred to as a virtual 
origin) is established for the turbulent analysis. 

For the hydraulically smooth case 
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Cf from equation (21) 
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For the distributed rough case 
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The turbulent flat plate method used here is that of Schlicting 
C8, chapter XXI) which is based on a transposition of Nikuradse*s densely- 
packed sand grain roughened pipe data* The effect of coiT5)ressibili-ty is 
due to the reduction in density at the wall as proposed by GoddardP 
The selection of the equivalent sand grain roughness for a given manufac- 
turing surface finish is made with the aid of Table II which was taken from 
Clutter! 2 
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TABLE II 


Equivalent Sand Roughness 
Type of Surface Kg (inches) 

Aerodynamically smooth 0 

Polished metal or wood 0.02 - 0.08 x 10"^ 

Natural sheet metal 0.16 x 10’^ 

Smooth matte paint, carefully applied 

Standard camouflage paint, average application 

Camouflage paint, mass -production spray 

Dip-galvanized metal surface 

Natural surface of cast iron 10 x 10"^ 


0.25 X 10"^ 
0.40 X 10'^ 
1.20 X 10'^ 

6 X 10'2 


Thickness Corrections 

The foregoing evaluations produce as estimate of the shearing forces 
on a flat plate (at zero angle of attack) for a variety of conditions. As 
an actual aircraft has a non vanishing thickness , an estimate of pressure 
gradient effects on skin friction and boundary layer displacanent pressure 
drag losses is required. A common procedure for acconplishing this and the 
one which will be used here is based on non- lifting experimental correlations 
for symmetric two-dimensional airfoils and .axisymmetric bodies. The 
following relations derived by Homer ( 13 , chapter VI) are used respectively. 


K 




1 K, -f ♦ 
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Homer recommends K, - 2 for airfoils with maximum thickness at 30% 
chord and K|= 1.2 for NACA 64 and 65 series airfoils. In this regard, the 
best information available to an analyst for his particular contour should 
be used. This is especially true for modem high performance shapes such 
as the supercritical airfoil, etc. 

Total Viscous Drag 

The aircraft total viscous drag coefficient is estimated by a sum of 
the preceding analysis over all components (i.e. wing, fuselage, vertical 
tail , etc . ) . That is 


M 

Co = 

viseooi 

4\l 





The con^jonent length used in the calculation of the skin friction 
coefficient is the mean chord for planar component segments and the physical 
length for bodies and nacelles. 


Zero Suction Drag 


The zero suction drag due to lift is calculated by numerically 
integrating the net pressure distribution times the projected area in the 
streamwise direction over each of the planer surfaces. The following 
block integration scheme is used to sura over all quadrilateral panels. 

c. • f; 2. Cp A. <x.; 

where 


and 
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cKq- is due to twist and camber, S is the control surface deflection 
and <ST= 1 for control surface panels and 0 for non control surface 
panels. Fj « 2 for symmetric geometries and Fi =» 1 for asymmetric 
geometries. For configurations having jet flaps, an additional drag term 
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must be included to account for the loss in thrust due to the downward 
deflection of the jet. This drag increment can be seen from the following 
thrust diagram. 



The section thrust loss can be expressed by the equation 

where 5^,. is the total deflection angle at the section of interest . The 
total increment can be found by a spanwise summation of the section thrust 
loss: 


AC 

D 



Potential Form Drag 


One hundred percent suction drag due to lift and st^ersonic wave drag 
due to thickness can be evaluated by integration of the momentum flux 
through a large circular cylinder centered on the x axis and whose radius 
approaches infinity. 


£ 
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The resulting expression for the total pressure drag is as follows: 


A, A, 

The first term represents the wave drag due to momentum losses thru 
the side of the cylinder caused by standing pressure waves. The second 
tenii represents the vortex drag which arises fran the kinetic energy left 
behind in the Treffiz plane by the system of trailing vorticies. Included 
in the vortex drag is the loss in the axial conponent of thrust due to a 
deflected jet flap. The calculation of these terms is discussed below. 

Vortex Drag 


The vortex drag may be conputed when the distribution of trailing 
vorticity in the Trefftz plane is known. The assunptions of linearized thin 
wing theory result in a vortex sheet which extends directly downstream of 
all lifting surfaces. By changing a surface integral for kinetic energy to 
a line integral over the vortex sheet in the Trefftz plane the following 
integral for drag results ' 
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where in addition 
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and where 


is the sectional 
is the sectional 

is the local normal velocity on the vortex sheet in the 
Trefftz plane 

9 is the local inclination of the sheet with respect to the 
y axis 
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The program ccanputes the normal velocity on the vortex sheet, , 

by assuming the vortex sheet is composed of finite trailing horseshoe 
vorticies whose strength is proportional to the local section C^(s). The 
normal velocity is computed at a control point located midway between the 
trailing vortex segments. 
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Each vortex j induces a contribution to the normal velocity at 
section i. 
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Wave Drag 


The integral for wave drag 




^CP 



r dlx <k9 


may be simplified by allowing the cylindrical surface of integration to 
recede infinitely far from the disturbance . Under these conditions 
it is possible to reduce spatial singularity simulations to a series of one- 
dimensional distributions. The basis for this reduction is the finding by 
Hayes (1^) that the potential and the gradients of interest induced by a 
singularity along an arbitrary trace on a distant control surface, say PP' 
of figure 6 (or alternately described by the cylindrical angle®), is 
invariant to a finite translation along the surface of a hyperboloid 
onanating from the trace and passing through the singularity. As the apex 
of the hyperboloid is a great distance away, the aforonentioned movement is 
along a surface which is essentially plane; it will be henceforth referred 
to as an "oblique plane". Since a singularity is a solution of a linear 
differential equation, all singular solutions which lie on the surface of the 
same hyperboloid (oblique plane) may thus be grouped to form a single equi- 
valent point singularity whose strength is equal to the algebraic sum of the 
individual strengths and which induces the same potential (momentum) along 
the trace as the group of individual singularities. 

This finding provides the basic technique for reducing a general spatial 
distribution of singularities to a series of equivalent lineal distributions. 
This is acconplished by surveying the three-dimensional distribution longi- 
tudinally at a series of fixed cylindrical angles, © . At each angle, the 
survey produces an equivalent lineal distribution by systematically cutting 
the spatial distribution at a series of longitudinal stations along its 
length. At each cut, the group of intercepted singularities is collapsed 
along the "oblique plane" to form one of the equivalent point singularities 
ccHEprising the lineal distribution . 
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The far field expression, for the wave drag of a general system 
of lift> and side force elements is 
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where 


■^0 3^.9) - 

;s'j = 


+(6,S) - 3 (S,©)^e - 3 te.s) c^e 

is the equivalent lineal singularity strength at the 
cylindrical angle 8 

equivalent source strength per unit length 

equivalent lifting el^ient strength per unit length 
equivalent side force strength per unit length 


These strengths are deduced from the three dimensional singularity 
distributions by application of the superposition principle along 
®^.^ipotential surfaces • For a distant observer such suirf aces are planar in 
the vicinity of the singularity configuration. The individual singularity 
strengths are related to the object under consideration by the requironent 
tangency at the solid boundary. Lomax Cl 5) derived the following 
approximate expressions between the equivalent singularity strengths and a 
slender lifting object. 
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where (see figure 7 

) 



A C«,^) is the Y-Z projection of the obliquely cut 

crossectional area 

c is the contour around the surface in the oblique cut 
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Figure 7 Areas and Forces Pertinent to the Evaluation of 
Wave Drag from the Far Field Point of View 




utilizing the singularity strength expressions derived by Lomax, the 
following expression for wave resistance based on the far field theory of 
Hayes is obtained ^ 
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In order to facilitate subsequent discussion, the above result is 
manipulated into the following form 
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where 




A(e,a) 


0 C 





A requirement for this transformation is that 
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In accordance with equation 23, the wave drag of a configuration is 
the average of the wave drag of a series of equivalent bodies of revoluation. 
The drag of each of these bodies is calculated from a knowledge of its 
longitudinal distribution of normal cross sectional area. For each 
equivalent body, these areas are defined to be the frontal projection of the 
areas and the accumulation of pressure force in the theta direction 
intercepted on the original configuration by a system of parallel oblique 
planes each inclined at the given Mach angle. The common trace angle ((|)) 
of the system identifies the equivalent body under consideration. 

Nacelles are assumed to swallow air supersonically. That is, the duct 
is operating at a mass flow ratio of unity. Consistent with this assun^tion, 
the equivalent body cross sectional area distribution is increased by the 
oblique projected duct capture area at all stations ahead of the duct 
which are intercepted by an oblique plane. 

Blunt base components are extended (maintaining constant cross 
sectional area) sufficiently far downstream to prevent flow closure around 
the base. 

In addition to a geometric description, a definition of the pressure 
distribution acting on the configuration is required. The vortex panel 
analysis is used for this purpose. The thickness pressures for planar 
components have tacitly been neglected under the assunption that the surfaces 
are sufficiently thin that the net pressure coefficient is representative of 
pressure acting on the oblique section. 

Estimation of the wave drag based on equation 23 depends on solution 
of integrals of the type 


I I 

O o 


of a numerically given function G(X). Evaluation of such forms has been 
studied by Eminton ^ ^ ^ 7 for functions having G ’ (X) continuous on the 
interval (0,1) and G’ (o) * (1) » 0. In such situations G^(X) can be 

expanded in a Fourier sine series. It can then be shown that 
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where 
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Eminton then solved for the value of the Fourier coefficients which result 
in I being a minimum, subject to the condition that the resulting series 
for GCX) be exact for an arbitrarily specified set of points o,l, ^i, i “ 1,^, 

This approach produces the following result 
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The solution of equation 23 for wave drag is accon^lished by use of the 
following identities. 
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Trim Drag 


The control surface deflection required to longitudinally balance the 
configuration at a given angle of attack is calculated from the pitching 
moment equation: 


Cr», Cm **■ ^TH1N\ =r O 
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(24) 
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The associated lift coefficient is given by 


« Cu iS- oc + I Cu J*. 


^ •< 


u 




(25) 


The trim drag is evaluated using the preceding drag analysis with the 
control surface deflection angle determined from equation 24. 

The lift and pitching moment characteristics appearing in equations 
(24) and (25) are based on the constant pressure vortex panel simulation 
described previously. Since this analysis is potential, dynamic pressure 
losses associated with viscous wake effects have been neglected. 
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CONCLUSIONS 


An aerodynamic configuration evaluation program has been developed and 
inqjlemented on a time sharing system with an interactive graphics teminal 
to maximize responsiveness to the preliminary analysis problem 

The solution is based on linearized potential theory and treats thick- 
ness and attitude problems at subsonic and st?)ersonic speeds. Three 
dimensional configurations having multiple non-planar surfaces of arbitrary 
planform and open or closed slender bodies of non- circular contour may be 
analyzed. Longitudinal and lateral -directional static and rotary derivatives 
solutions may be generated. 

Nominal case computation time of 45 CPU seconds on the CDC 175 for 
a 200 panel simulation indicates the program provides an efficient analysis 
for systematically performing various aerodynamic configuration tradeoff 
and evaluation studies. 
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